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Summary. We prove a hydro dynamic limit for ballistic deposition on a multidi- 
mensional lattice. In this growth model particles rain down at random and stick 
to the growing cluster at the first point of contact. The theorem is that if the 
Q\ . initial random interface converges to a deterministic macroscopic function, then at 

^ ! later times the height of the scaled interface converges to the viscosity solution of 

1-^ I a Hamilton-Jacobi equation. The proof idea is to decompose the interface into the 

^ • shapes that grow from individual seeds of the initial interface. This decomposition 

converges to a variational formula that defines viscosity solutions of the macrosopic 
equation. The technical side of the proof involves subadditive methods and large 
K> ', deviation bounds for related first-passage percolation processes. 
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1. Introduction 

In a ballistic deposition model particles fall on a surface, find a location where 
they attach, and become part of the growing cluster. Depending on the rules chosen, 
the particle may stick to the first point of contact, or it may "roll downhill" and 
attach itself to the first stable location it encounters. We study the version where 
particles stick to the first point of contact. A ballistic deposition model is flux 
limited as opposed to reaction limited because the rate of growth is limited by the 
availability of material rather than by the availability of growth sites. 

Another qualitative distinction between growth models is the division in /oca/ and 
nonlocal models. An example of a local model is one where the rate of attaching 
new particles depends only on the states of some finite number of neighboring sites. 
A highly nonlocal model is diffusion limited aggregation where the chances of a 
particle attaching at any particular site are influenced by far-away parts of the 
growing cluster. Ballistic deposition has a nonlocal aspect because it can build 
overhangs that extend far sideways to shade parts of the interface that then no 
longer receive the particle flux. However, the height of the interface has local 
dynamics, and it is this height process that we study in our paper. 

The shadowing effect and low atomic mobility (stick to the flrst point of contact) 
may play a role in actual deposition processes, so there is physical motivation for 
these rules (see [16]). We refer the reader to the survey article of Krug and Spohn 
[18] for a general discussion of growth models and for references to the physics 
literature. 

We prove a hydrodynamic scaling limit for a ballistic deposition process on an 
inflnite cZ-dimensional cubic lattice. Particle deposition happens independently at 
all sites, governed by exponential waiting times. Once deposited, particles never 
leave the surface, so in particle system jargon this process is totally asymmetric. The 
correct hydrodynamic scaling shrinks space and speeds up time by the same factor 
n. The theorem is a law of large numbers: as n -^ oo, the height of the randomly 
evolving surface converges to the solution of a flrst-order partial differential equation 
of the Hamilton-Jacobi type. For general treatments of hydrodynamic limits, the 
reader may consult references [5], [15] and [29]. 

From a statistical mechanics point of view, our result is a rigorous microscopic 
derivation of the expected macroscopic theory. The macroscopic growth velocity 
of the height function -0 is determined by the local slope: dt/j/dt = /(V^). The 
function / that gives the dependence is the Legendre conjugate of the stationary 
shape (7 of a cluster grown from a seed: f{p) = sup^{x ■ p + g{x)}. As explained in 
[18], this can be viewed as a Wulff construction for a growing shape. 

The essence of the proof is to construct the process so that a supremum of 
ballistic deposition processes is again a ballistic deposition process. The proof 
works equally well for discrete time in which case the waiting times of deposited 
particles are geometric rather than exponential. In fact the Markovian nature of the 



dynamics is not needed at all. The proof can be carried out for arbitrary waiting 
times, but we have not bothered with such generality here. An example can be 
found in [26]. 

The paper is organized as follows. Section 2 defines the model more precisely and 
states the theorems. We discuss briefiy viscosity solutions of Hamilton- Jacobi equa- 
tions, and suggest open problems from the physics literature that may be amenable 
to rigorous progress. In Section 3 we do the standard graphical construction (see 
for example [6], [9], [11], or [19]) to show that the process can be rigorously defined 
on a probability space that supports the Poisson jump time processes. And we 
construct the coupling that expresses the process as the supremum of the shapes 
from individual seeds of the initial interface. 

In Section 4 we prove exponential large deviation estimates for a first-passage 
site percolation process. We need these because the lateral (sideways) growth of 
ballistic deposition from a seed is equivalent to a first-passage percolation problem. 
The bounds we need are of a standard type, but they are not available in the 
literature in the right form. We use a result of Talagrand [30] and some ideas from 
Kesten [14] to derive the inequalities. Grimmett and Kesten [10] have proved such 
bounds for growth along coordinate axes, but we need to control growth in all space 
directions. Generalizing the block argument of [10] appeared harder than applying 
the very general tools of [30]. 

Theorem 1 is proved in Section 5, and in Section 6 a technical extension of 
Theorem 1 is proved where the initial seed and initial time are translated as the 
limit is taken. This is required for the proof of Theorem 2, which is the content of 
Section 7. 

Notational conventions. Frequently used notation is summarized here for the 
reader's convenience. C, Ci, C2 are constants whose exact values are immaterial 
and whose values may change from one inequality to the next. R-|- = {x G R : 
x > 0} is the set of nonnegative real numbers, and similarly for Q_|_ and Z+. 
N = {1, 2, 3, ... } is the set of natural numbers. Sites or points of Z'^ are denoted 
by w, Vj w, z. is the origin of Z'^, and J\f is the set of nearest neighbors of the 
origin in Z'^, in other words, the 2d sites at £^-distance 1 from the origin. To 
distinguish points of Z'^^^ from those of Z'^, we call points of Z'^^^ cells and denote 
them by (w, k) where w G Z'^ and k E Z. In particular, (0, 0) is the origin of 
Z'^ X Z_|_. For a real number x, [x] is the maximal integer n subject to n < x. 
For X = (xi,... ,Xd) G R'^ the site [x] = ([xi],... , [x^]), and the i^ norm is 
\x\ = |xi| + ■ ■ • -|- \xd\- For X, y G R, X V y = max{x, y}, x A y = min{x, y}. 
\\g\\oo = SUP3, \g{x)\ for any function g. 

1a and 1{A} denote the indicator random variable of the event A. An Exp(/?) 
random variable X satisfies P(X > t) = exp(— /3t) for t > 0. S^ stands for a sum 
of n i.i.d. Exp(l) random variables. The standard notion of stochastic dominance 
is expressed hy X ^ Y which is equivalent to P{X < t) > P{Y < t) for all t. Og 



is a time translation on Poisson point processes that translates points r to r — s. 
In other words, reading OgUJ from time onwards is the same as reading u from 
time s onwards. Bq C R'^ is the convex compact limiting shape for first-passage 
site percolation on Z*^ with Exp(l) waiting times. 



2. The results 

We start with a description of the ballistic deposition process. Fix a positive 
integer d. Imagine first an arbitrary subset of Z*^"*"^ of occupied cells. Cells are 
simply points of Z'^^^ and denoted by (w, h), where -u G Z*^ is a site and /i G Z is a 
height. The cluster (subset of occupied cells) grows through deposition events: over 
each site tt G Z'^, particles rain down randomly at exponential rate 1, independently 
of all other sites. A particle descending down sticks to the first spot where it touches 
the existing cluster. We imagine that the particles are exactly the size of a unit 
cube in d + 1 dimensions. Over a site u a particle instantaneously drops down from 
height /i = oo to the highest cell adjacent to the existing cluster, and then this cell 
becomes occupied and joins the cluster. (Cells are adjacent if their i^ distance is 
1.) As the process evolves, a porous structure is generated that grows upward in 
d+1 dimensions. 

We follow the evolution of the top surface of the cluster. For each site u, the 
height au is the maximal h such that cell (-u, h) is occupied. We permit the values 
au = ±oo. au = — oo means that no cell in the column {{u,h) : h E Z} is occupied. 
The state of the process is the configuration a = {au ■ u G 7/^) of height variables, 
and the state space is 5 = (Z U {±oo})'^ . We write a{t) = {au{t)) for the height 
process, where t > denotes time. 

The time evolution of the interface is determined by a collection T = {T" : u G 
2^1 Qf y-Q^Q I Poisson point processes on the time line (0, cxo). At each epoch of T" a 
particle is deposited above site u. In terms of the height variables au{t) a deposition 
event has a simple expression: Let A/" denote the set of 2d nearest neighbors of the 
origin in Z*^. If r is an epoch in Poisson process T", the height above site u jumps 
at time r according to this formula: 

(2.1) au{r) = max[au{r-) + 1 , ma.x{au+z{r-) : z G A/"}] . 

The maximum in (2.1) has the effect that the deposited particle sticks to the highest 
unoccupied cell above site u that is adjacent to the existing cluster. If no cell above 
site u is adjacent to the cluster at time r— , the right-hand side of (2.1) equals — oo. 
This means that the deposited particle is lost because it cannot stick to the cluster. 
In Section 3 we construct this process rigorously starting from an arbitrary initial 
interface cr(0). 

The simplest ballistic deposition process starts from a single occupied cell, or 
seed. The prototypical one starts from a seed at the origin. For this process we 



reserve the symbol Z and use the symbol a for the general ballistic deposition 
process. So at time 0, ^o(O) = 0, and ^u(O) = — oo for all sites -u 7^ 0. Otherwise 
Z'(-) evolves as specified above. 

Let i3(t) denote the set of sites above which Z has an occupied cell by time t: 

(2.2) E{t) = {w G Z^ : Z^it) > 0}. 

Since the original seed is at the origin, Zu{t) > is equivalent to Zu{t) > —00. The 
cluster B{t) is the view of the process Z{t) from above, by projecting Z^"*"^ onto 
Z*^ by the map (ti, h) 1-^ u. The ballistic deposition rules imply that B{t) grows 
according to this rule: Each site adjacent to the existing cluster joins independently 
at rate 1. [Because particles rain down at rate 1 above each site u, and they have 
a chance of sticking iff some cell above an adjacent site is already occupied.] Thus 
B{t) is a familiar growth model, namely first-passage percolation, or a version of the 
Eden growth model. It satisfies a law of large numbers. There is a closed, convex 
deterministic set Bq C R^ with nonempty interior such that this holds almost 
surely: for any £ > 0, 

(2.3) [t{l - e)Bo] n Z'^ C B{t) C [t(l + e)Bo] n Z'^ 

for all large enough t. In first-passage percolation literature the random waiting 
times are usually attached to the edges, but for our B{t) the waiting times are 
attached to the sites. In Section 4 we prove some large deviation estimates for 
B{t). 

Over the set Bq the scaled ballistic deposition process t ^Z{t) approaches a 
limiting shape. Let intSo denote the topological interior of Bq. We prove in 
Section 5 this theorem: 

Theorem 1. There exists a bounded positive function g defined on the open convex 
set intSo such that this law of large numbers holds: Outside an event of probability 



(2.4) lim -Z[^,^{nt)=tg{x/t) 

for all X G R'^ and t > such that x/t G intSo- Furthermore, g is continuous, 
concave, and invariant under permutations of the coordinate axes and reflections 
about the origin. 

The main theorem is a scaling limit from a general initial interface. Suppose 
we have a sequence cr" of ballistic deposition processes. The initial macroscopic 
interface is given by a function -i/'o on R'^. We consider three different sets of 
hypotheses. 



Assumption A. i/jq is a continuous [—00, +oo]-valued function on R*^. The ballis- 
tic deposition processes o""(t) are constructed on a common probability space, and 
all processes use the same version of the Poisson jump time processes. There is a 
countable dense subset Yq C R'^ such that 

(2.5) lim la[;^](0)=^o(?/) 

almost surely for each y & Yq, and each y ^ Yq has arbitrarily small closed neigh- 
borhoods V such that almost surely 

(2.6) lim sup — ■ max o"^(0) < sup -00 (y) . 



Assumption B. i/^o is an arbitrary [—00, +oo]-valued function on R'^. The ballistic 
deposition processes cr'^(t) are constructed on a common probability space, and all 
processes use the same version of the Poisson jump time processes. For any fixed 
y G R'^, the limit (2.5) holds almost surely, and there are arbitrarily small closed 
neighborhoods V oi y such that (2.6) holds almost surely. 

Assumption C. Again ■i/'o is an arbitrary [— cxd, +oo]-valued function on R^. (2.5) 
and (2.6) hold in probability in this sense: Let (O^, JF„, P^) be the probability space 
of the nth process cr". For any y E H'^ and £ > 0, 

(2.7) lim P^(\n-^a\'^y^{0)-My)\>e) =0, 

and there are arbitrarily small closed neighborhoods V of y such that 



(2.8) lim Pn max o-"(0) > n ■ sup ipo{y) + ne ] = , 

n— >oo \ uEnV yEV 



For X G R and t > 0, define 



(2.9) ^{x,t)= sup <'^o(l/)+t(7(^^ 



and V'(a^7 0) = '^q{x). 



Theorem 2. (a) Strong Laws of Large Numbers: Under Assumption A we have this 
convergence, simultaneously for all x G R'^ and t > 0, outside a single exceptional 
event of probability zero: 

(2.10) lim ia[;,](nt)=^(x,t). 

Under Assumption B the limit in (2.10) holds almost surely for all (x^t) at which 
ifj is t-continuous from the right, ^^{x, t) = '^{x, t+). 

(b) Weak Law of Large Numbers: Under Assumption C the limit in (2.10) holds 
in probability for all (x, t) at which ifj is t-continuous from the right. In other words, 
for all such (x, t), 

(2.11) lim P^(\n-^aV.{nt)-i;{x,t)\>e] =0 

n — >oo \ I 

for all e > 0. 



Remarks and extensions 

2.1 The assumptions. Some uniformity assumption such as (2.6) is needed for the 
resuh. Consider this example in dimension d = 1: The initial interface is 

f n, u = 1 
"^ ^ I 0, w^l. 

The limit (2.5) is satisfied with il^o{x) = 0, and then from (2.9) 'i/j{xjt) = tg{0). But 
the spike in a""(0) at site u = 1 implies that 



n 



~^ao{nt) > 1 + n ^ZQ{nt) — > 1 + tg{0) as n ^ oo, 



where we wrote Z^ for the process that starts from a seed in cell (1,0). 

The point of separating Assumptions A and B is that continuity of i/jq guarantees 
continuity of '0(x, t) [Lemma 7.1 in Section 7]. Then the limit (2.10) can be proved 
almost surely, simultaneously for all (x, t), even though the hypothesis requires the 
limit only on a countable dense set Yq. The limit under Assumption B generalizes 
Theorem 1. The right-continuity in t at {x,t), namely il){x,t) = 'i(^{x,t+), follows 
from x-continuity at (x,t), as can be verified from (2.15) below. 

The nonlocal shadowing effect restricts the law of large numbers to continuity 
points of the macroscopic surface. Consider the process Z growing from a seed at 
(0,0) in dimension d = 1. The limit (2.4) cannot possibly hold at a; = t because 
whether site [nt] is shadowed by the top surface varies with the fiuctuations of the 



right edge of the growing cluster: Z'[^i](?it) = — oo or > depending on whether 
'^int] > nt or < nt, where SI^t^ is a sum of [nt] i.i.d. Exp(l) random variables and 
represents the time it takes for the cluster to grow [nt] lattice units to the right. 

The proof of Theorem 2 depends on two essential things: a coupling (Lemma 3.3) 
that expresses the general process as a supremum of processes of the Z type, and 
the limit in Theorem 1 for the Z process. The coupling is basically a combinatorial 
property of the paths. As long as pathologies such as jump times accumulating at 
a point are ruled out, the coupling is independent of the probability distribution 
on the paths. The limit in Theorem 1 on the other hand comes from an ergodic 
theorem and is closely tied to the probabilities of the evolution of Z. Thus our proof 
of Theorem 2 works for any reasonable jump rates under which Theorem 1 can be 
proved. For example, we could take quenched random jump rates where deposition 
at site u happens at rate Uu, and the rates {a^} are i.i.d. random variables (for 
examples in other contexts, see [2], [26] and [27]). Another possibility would be to 
let a sufficiently regular function a{x) on R'^ determine the rates so that deposition 
at site u for process a" happens at rate a{u/n) (for examples, see [1] and [21]). 

2.2 Viscosity solutions of Hamilton- J acohi equations. Consider the Hamilton- Jacobi 
equation 

(2.12) ^t - fiy^) = 

(2.13) '4^{x,Q)=Mx) 

where / is some continuous function on R'^. A function 'i/j{x,t) on R*^ x R_|_ that 
satisfies the initial condition (2.13) is a viscosity solution of (2.12)-(2.13) if the 
following holds for all continuously differentiable functions (j) on R'^ x (0, oo): if 
tl) — cj) has a local maximum at (xq, to), then 

(/)t(xo,to) - /(V(^(xo,to)) < 0, 

and ii if) — (f) has a local minimum at (xq, to), then 

</'t(a;o,to)-/(V(/)(xo,to)) >0. 

The notion of viscosity solution is due to Crandall and Lions [4]. Properties of 
viscosity solutions of Hamilton- Jacobi equations are developed in Crandall, Evans 
and Lions [3] , and in the textbook of Evans [8] . 

Equation (2.9) is known as a Hopf-Lax formula in the partial differential equa- 
tions literature. Let / be the (negative of the) Legendre conjugate of g: 

(2.14) f{u)= sup {u-x + g{x)}, u e K"^ . 

xGint Bf) 



Since Bq is a compact set and g is bounded, one can check that / is finite and 
Lipschitz continuous on aU of R^. Then an argument in [8] (proof of Theorem 3 in 
Section 10.3.4) shows that ifj^x^t) defined by (2.9) is a viscosity solution of (2.12)- 
(2.13). By applying results from the p.d.e. literature, we can refine Theorem 1 with 
a uniqueness statement: 

Theorem 3. Suppose Assumption A is in force, and that additionally the initial 
macroscopic profile i/jq is uniformly continuous on R^. Fix a finite time horizon 
T < oo. As in Theorem 1 the strong law of large numbers (2.10) is valid. On 
R'^ X [0, T] the limit i/j{x,t) is the unique uniformly continuous viscosity solution 
of the Hamilton- Jacobi equation (2.12)-(2.13) whose Hamiltonian —f is defined by 
(2.14). 

Proof of Theorem 3. The point is that the additional assumption of uniform con- 
tinuity on ^/'o permits us to make the uniqueness assertion about '0. Uniqueness 
theorems for unbounded viscosity solutions have been proved by Ishii [12], whose 
Theorem 2.1 states that equation (2.12)-(2.13) with continuous / has a unique uni- 
formly continuous viscosity solution on R'^ x [0, T]. Assume Theorem 2. To prove 
Theorem 3, we therefore need to check the uniform continuity oiij) defined by (2.9), 
assuming that %1)q is uniformly continuous. Here is an outline of the argument. 

First check that formula (2.9) operates like a semigroup: Once tp has been defined 
by (2.9), it follows for all < s < t that 



(2.15) ^{.x,t)= sup \^{y,s) + {t-s)g['^-^ 

yex+{t-s)-intBo I \l — S 



This bounds the growth of i^{Xjt) in time. Let bo = sup{|y| : y G Bq}. Fix x and 
s < t. (2.15) gives 

< ■i/j{x,t) — i{j{xj s) 

(2-16) < sup \^^y'^s)-ij{y",s)\ + it-s)\\g\\^. 

\y'-y"\<bo(t-s) 

On the other hand, for a fixed time t the definition (2.9) directly bounds the vari- 
ation of -(/;(• , t): For any Xi, 0:2 G R^ and t > 0, 



(2.17) |^(a;i,t)-v(x2,t)| < sup \My') - My")\ 

\y'-y"\<\xi-x2\ 
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Now the uniform continuity of -0 on all of R^ x R_|_ follows: 

\ij{xi,ti) -tp{x2,t2)\ 
< ll/jixi^ti) - ij{x2,ti)\ + |^(x2,ti) - V'(a^2,t2)| 

< sup \i^oiy') - ipo{y")\ 

\y'-y"\<\xi-x2\ 

+ sup \i^{x',ti)-ij{x",ti)\ + {t2-ti)\\g\\oo 

\x'-x"\<bo(t2-ti) 

[ without loss of generality, assume that ti < ^2 ] 

< 2 ■sup{|V'o(y') -i^o{y")\ ■ \y' -y"\ < \xi -X2\ v6o|^2 -^i|} 

+ |t2- til -11(71100. 

This proves Theorem 3. D 

2.3 Statistical mechanics. From the point of view of statistical mechanics, our 
paper provides a rigorous derivation of the macroscopic theory that is taken as basic 
in the physics literature. According to this theory, macroscopically the interface 
moves under an inclination-dependent growth velocity f{u), and / is the Legendre 
conjugate of the cluster shape g{x) that grows from a seed. From this basis the 
physics literature seeks to describe finer properties of the deposition process. The 
reader is referred to the survey article [18], and to articles [16], [17], and [20]. 

Here we comment on some properties of the macroscopic objects, and mention 
open problems suggested by the physics papers. In general, describing / and g 
is as hard as first-passage percolation, since knowing g would imply knowing the 
first-passage percolation shape Bq. In dimension d = 1 the percolation question 
becomes trivial, So = [—1, 1], so one may hope to say something about / and g for 
d = 1. The case rf = 1 is also the focus of the physics literature. 

It follows from (2.14) that the velocity f{u) is convex and even [f{u) = f{—u)]. 
Consequently it has a minimum at w = 0, /(O) = g{0). Whether this minimum is 
strict as expected (p. 903 in [17]) is a harder question because that would require 
showing that g does not have a corner at a; = 0. These types of questions are 
open for other interface models too, except in those rare cases where invariant 
distributions can be used to explicitly compute limiting shapes. For examples, see 
[23]~[27]. In d = 1 (2.14) gives linear asymptotics f{u) = \u\ + g{l) + o{\u\) for 
large slopes u — > ±oo. The precise nature of the error 0(1^1) would be of interest. 
Equivalently, one wants to know the asymptotics of g{x) — g{l—) as x /" 1. 

The velocity must increase with dimension, since higher dimension means more 
neighbors to speed up the growth over a particular site. This is easy to check by 
coupling the Z-processes for dimensions d and d + 1 so that the (i-dimensional Z- 
process grows over the hyperplane {xd+i = 0} in Z*^"*"^. Without any estimates, this 



11 
gives gd+i{x'jO) > gd{x') for any x' G R'^. Consequently, for any u = {u',Ud+i) G 

fd+i (u) > sup {u ■ X + 5fd+i (x', 0) } 

> sup {a -x ^ gd{x'y\ 

x'eR'' 

= h{y!) ■ 

Mean-field analysis in [17] suggests that /d(0) grows like (i/log(i, and simulations 
appear to show a slow convergence toward mean- field values as (i ^ oo. Growth at 
rate (i/log(i has been verified for first-passage percolation [13], so these questions 
can be investigated rigorously. 

In a synchronously updated ballistic deposition process time is discrete, and a 
rate parameter p G (0, 1) is fixed. At each time step t = 1, 2, 3, . . . , an independent 
random choice is made at each site: with probability p the height is updated ac- 
cording to equation (2.1), and with probability 1 — p the height remains the same. 
The results of our paper apply to the synchronous process without changes. The 
only difference is that the Poisson point processes {T™ : w G 7/^} of jump times 
are replaced by Bernoulli processes on the discrete time line N = {1,2,3,...}. In 
these processes an event arrives at each time with probability p, independently of 
everything else. Analogously with the flat edge result of Durrett and Liggett [7] 
for first-passage percolation, a faceting transition happens in ballistic deposition for 
large enough p. Interestingly, simulations in [17] suggest that the velocity f{u) is 
linear for all slopes -u > 1 if p is large enough. 



3. Construction and coupling 

To construct the ballistic deposition process, start by giving each site u E Z'^ a,n 
independent rate 1 Poisson point process T" on the time line (0, oo). Fix an initial 
configuration a(0) = {(Ju{0) : u G Z*^) G iS = (Z U {±oo})^ . Informally speaking, 
the construction of the dynamics goes as follows: If r is an epoch (in other words, 
a point) of T", then at time r height variable a^ jumps: 



(3.1) au{r) = max[au{r-) + 1 , max{au+z{r-) : z G A/"}] . 

Recall that J\f is the set of nearest neighbors of the origin in Z'^. 

To make the construction rigorous, we show that there exists a fixed time to > 
and a set of Poisson processes {T^} of full probability such that, starting with an 
arbitrary a{0) G 5, the evolution a{t) can be computed for t G [0,to]. Since to 
is independent of a"(0), the construction can be repeated, starting with cr(to), and 
extended to time interval [0, 2to] . And so on, to arbitrarily large times. 
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Given a fixed number to > and the Poisson processes {T"}, construct the 
foUowing random graph with vertex set Z^: Connect nearest neighbors u and v 
with an edge if either T" or T" has an epoch in [0, to]. 

Lemma 3.1. For small enough fixed to > 0, this random graph has no infinite 
connected components for almost every realization of {T"^}. 

Before proving the lemma, let us see how the construction problem is solved. 
Make these further assumptions on the Poisson point processes, valid for almost 
every realization: 

(3.2a) The {T"} are such that there are no simultaneous jump attempts. 
(3.2b) Each T^ has finitely may epochs in each bounded time interval. 

All sites w that can infiuence the evolution at site u up to time to are connected 
to u in the random graph. Since u lies in a finite connected component C, the 
point process UwecT^ has only finitely many epochs in [0,to]. Consequently the 
evolution au,{t) can be computed for zw G C and t G [0, to] from rule (3.1), by 
considering the finitely many epochs in their temporal order. This procedure is 
repeated for all connected components. 

Proof of Lemma 3.1. By translation invariance, it suffices to show that the origin 
is almost surely connected to only finitely many vertices, {uq, -Ui, . . . , -u^} is a self- 
avoiding path of length n in the random graph if Ui ^ Uj for i ^ j and if there 
is an edge between Ui and Ui+i for each i. If the origin is connected to a site u 
with 1^1 > L, there is a self-avoiding path of length > L starting at the origin. The 
probability that a self-avoiding path of length 2?i — 1 starts at the origin is at most 

(2rf)2-i(l-e-2*")". 

The factor (2(i)^"~^ is an upper bound on the number of such paths. If = 
uo,ui,. .. ,U2n-i is such a path, the n edges {uo,ui), (^2,^3), . .., (w2n-2, W2n-i) 
are present independently of each other, and each with probability 1 — e~^*" [at least 
one of T^2^ and T^^^+i must have an epoch in [0, to], and each T"^ has rate 1]. Pick 
to small enough so that (2(i)^(l — e~^*°) < 1. Then by Borel-Cantelli, self-avoiding 
paths from the origin have a finite upper bound on their length, almost surely. D 

This approach to the construction of a particle system is due to Harris [11]. Our 
presentation followed [6]. 

Let (O, JF, P) denote the probability space whose sample point to represents a 
realization of the Poisson processes T = {T^}. We constructed the random path 
^(■) — i'^ui't) '■ u £ Z^,t > 0) as a function of the given initial state cr(0) and a 
sample point u. Since the Poisson processes are Markovian, the process cr(-) is a 
time-homogeneous Markov process. When the initial interface (t(0) is random, the 
underlying probability space is constructed so that a{0) and T are independent. 

Formula (3.1) preserves ordering, so we get the following monotonicity lemma, 
whose proof is left to the reader: 
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Lemma 3.2. Suppose a and p are ballistic deposition processes constructed on a 
common probability space so that they use the same version {T^} of the Poisson 
processes. Assume that at time 0, C7u{Q) > P-u(O) for all u E Z'^, almost surely. 
Then almost surely au{t) > Pu{t) for all u eTa'^ and t > 0. 

A less obvious property of the construction is the foUowing, which forms the 
basis of our approach to the hydrodynamic hmit. 

Lemma 3.3. Suppose the ballistic deposition process a and a countable family of 
ballistic deposition processes {Q : z G X} are constructed on a common probability 
space so that they all use the same version {T^} of the Poisson processes. Assume 
that at time 0, almost surely, 

(3.3) a„(0) = supC(O) for all u E Z^. 

Then almost surely 

(3.4) auit) = supC(t) for all ueT^ andt>0. 



Proof. First apply Lemma 3.2 with p = C* to get 

(3.5) a^(t) > supC(t) for all w G Z'' and t > 0. 

iex 

Thus we need to show that for all sites u and times t, there is some index i such 
that(T„(t) =C(^)- 

Pick and fix a realization {T"} that satisfies assumptions (3.2) and for which 
the conclusion of Lemma 3.1 holds. We first show that, for any processes a and 
{Q : i E X} that satisfy the hypotheses, (3.4) holds for t E (0, to] where to > is 
the number chosen in Lemma 3.1. To do so for a fixed site w*^, let C C Z'^ be the 
finite connected component of u^ in the (random) graph constructed for Lemma 
3.1. The evolutions of all the processes on the sites of C are determined by the 
finitely many Poisson points in U^^,gcT^ fl (0, to]. We can now prove (3.4) up to 
time to by checking that it holds right after each jump. 

So suppose r G (0, to] is a jump time in T'" for some site v E C so that (3.1) 
happens for u = v. Assume by induction that (3.4) holds for t < r, for all u E C. 
Since the variables cr^ and Q^ are Z-valued, this means that at each time t < r the 
supremum in (3.4) is actually achieved at some i EX. Depending on how the jump 
(3.1) for -u = f is realized, two cases need to be considered. 
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Case 1. First suppose a^{r) = ay{r—) + 1. By induction, there is a j G X such 
that ay{r—) = C^(r— ). Since Q jumps too and by (3.5), 

Mr) >supC(r) >Ci{r) 

= niax[C^(r-) + 1 , max{C^_^^(r-) : z E Af}] 
>Q{r-) + l = ay{r-) + l 
= o-„(r). 

Thus Q jumps to the same height as a„. 

Case 2. The second possibihty is that for some w G A/", ay{r) = cr„_|_^„(r— ). Then 
by the jump rule (3.1) 

ay+w{r-) > ay{r-) + 1 and cr„+^(r-) > ay+z{r-) for all z E M. 

By the definition of the random graph v + w E C, so hy induction there exists a 
j eT such that ay+y,{r-) = C^+^(r-). By (3.5) 

(Jv{r-) + 1 > a(r-) + 1 and a,+,(r-) > C'+.(r-) for zEAf. 

Together these equalities and inequalities yield 

Ci+Jr-) > max[C^(r-) + 1 , max{C^+,(r-) : z E M}] , 

so the jump rule (3.1) applied to Q gives 

Q{r) = C+^(r-) = ay+y,{r-) = (7„(r). 

In other words, Q jumped to the same level as cr„, and (3.4) continues to hold right 
after the jump time r. 

The summarize: we have shown that (3.4) holds for times < t < to for any 
processes cr, {C} that satisfy the hypothesis (3.3) at t = 0. Now apply the same 
step again, to the processes a{t) = a{to+t) and (C^^) = C*(^o+^)}- These processes 
satisfy the hypothesis at t = by virtue of (3.4) at t = to. This way the validity of 
(3.4) is extended to times < t < 2to. And so on, to arbitrarily large times. D 

For f G Z"^, let Z"" = {Z^{t) : u E Z'^) denote the ballistic deposition process 
started from a seed in cell {v, 0). In other words, initially 

( if u = V 

(3.6) Z:{0) = \ 

I — OO 11 U ^ V . 
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Given an arbitrary initial configuration cr(0) = (cr^i(O) : u G Z'^) [random or deter- 
ministic], define a family of processes, indexed by X = Z'^, by the initial conditions 

w ^ ( O-yiO) if U = V 

(3.7) C(0) = ., . 

[ — OO 11 U ^ V . 

We can write 

(3.8) CJt) = a^iO) + Z:{t) 

with the convention that oo + (— oo) = — oo. The processes a and {C" : v E Z'^} 
satisfy (3.3). Lemma 3.3 gives this corollary, which is basic for our proof of the 
hydro dynamic limit: 

Corollary 3.1. The equality 

(3.9) au{t)= supK(0) + Z^(t)} 

holds almost surely, for all u E Z'^ and t > 0. 



4. First- PASSAGE site percolation 

In this section we prove some estimates for the first-passage percolation problem 
briefiy encountered in the introduction. First we redefine it in the standard way. 

Give each site -u G Z'^ an Exp(l)-distributed random time t('u), independently of 
the other sites. Say tt = {w^, w^, . . . , w"^} C Z'^ is a nearest-neighbor path from u to 
V of length m if to < oo, w° = -u, w"^ = v, and |w* — ty*~^| = 1 for i = 1, . . . ,rM. We 
use I ■ I to denote the l^ norm: \w\ = \wi\ + ■ • ■ + \wd\ for w = {wi, W2, . . . , Wd) G Z'^. 
The passage time of a path tv = {w^, w^, . . . , w"^} is 



(4.1) T{n) = J2t{ 



m 



Since the path starts from site w^, the value t{w^) is not included in the sum. The 
passage time from site -u to f is 

(4.2) T{u,v)=MT{n) 

TT 

where the infimum ranges over nearest-neighbor paths tt from u to v. The mini- 
mization has the effect that n may be assumed self-avoiding in the sense that there 
are no repetitions among {w^, w^, . . . , w"^}. 
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The cluster growing from a seed at the origin is defined by 

(4.3) B{t) = {ueZ'^ :T{0,u)<t}. 

It is clear from the description that B{0) = {0}, and B{-) grows according to this 
local rule: Each site adjacent to the current cluster joins independently with rate 
1. 

To make the connection with ballistic deposition, consider again the process Z 
started from a seed at the origin: 

.X . X r 0, w = 

(4.4) ^.(0) = 

[ — oo, u ^ . 

Let R{u, h) denote the first time Z is at or above height h E Z_|_ over site u G Z^: 

(4.5) R{u, h) = inf{t > : Z^{t) > h}. 

In particular, R{u, 0) is the first time that a particle sticks to the cluster above site 
u. -R(0, 0) = by definition (4.4). Since the notation may lead to confusion, let 
us emphasize that the of R{u, 0) is the zero of Z_(_, while the of T(0, u) is the 
origin of Z^. 

Lemma 4.1. We have the following equalities in distribution between the cluster 
and passage time processes: 

(4.6) l{ueZ'^:Zuit)>0}:t>o\ = {B{t) : t > O} 

and 

(4.7) {R{u, 0):ueZ'^} = {T(0, u) : u e Z'^} . 

Proof. Abbreviate B{t) = {u e Z'^ : Z^it) > 0} for this proof. B{-) and B{-) 
are both Markov jump processes on the countable state space of finite, connected 
subsets of Z'^ that contain 0. They have the same initial state B{0) = B{0) = {0}, 
and both processes have identical infinitesimal rates: each site adjacent to the 
current cluster joins independently at rate 1. Hence the two processes are equal in 
distribution, which is (4.6). 

By definition (4.3), T(0, -u) is the first time when site u joins the cluster B{t). 
R{u, 0) has the same meaning for B{t). So the processes {T(0, u)} and {R{u, 0} are 
obtained by applying a certain measurable function to the processes B{-) and B{-). 
Thus (4.7) follows from (4.6). D 
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In later sections (4.6) will be used several times to derive deviation bounds for 
Z. 

Return to T{u^v) as constructed by (4.3). Subadditivity considerations give a 
limit 

(4.8) lim — T(0, \nx\) = fx{x) a.s. 

for all X G R^. The function fx is convex, homogeneous [ii{rx) = rn{x) for r > 0], 
and Lipschitz continuous. The limiting cluster is defined by 

(4.9) Bo^{x eK'^ : ^i{x) < 1} , 

and now the set convergence (2.3) is valid. We shall not give proofs of these laws of 
large numbers. The case where the random times t{u) are on the edges instead of 
on the sites is thoroughly discussed in Kesten's lectures [13]. We prove some large 
deviation estimates for T(0, u) that we need in the sequel. 

Proposition 4.1. For any x G R'^ and £ > there are finite constants Ci = 
Ci{x,e) > such that 



P(|T(0, [nx])-n^{x)\ > ne) < Ci exp(-C2n) 

for all n. 

Before the proof we derive two corollaries, assuming Proposition 4.1. 

Corollary 4.1. For any t, e > there are finite constants Ci = Ci{t,e) > such 
that 

J2 P{T{0, u) < nt) < Ci exp(-C2n) 

weZ'* ,u(^n{t+e)Bo 

for all n. 

Proof of Corollary 4-i- Fix a large n and small Sq, Si, £i > so that 

(4.10) 6o + d/n < 6i and 6i{l + £i) < £ - £i. 

Let A^ = [{m + l)(t + e)Bo \ m{t + e)Bo] n Z^ for m > n. Write 

oo 

(4.11) Yl PiTiO,u) <nt)=J2 E ^(^(0'^) ^ ^0 • 

u^n{t+e)Bo m=nueAm 

Since \Am\ < Cm'^, it suffices to bound P{T(Q,u) < nt) < exp(— Cm) for u G Am- 
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Pick and fix points x^, . . . , x'^ G 2(t + e)Bo \ {t + e)Bo such that each y G 
2{t + e)Bo\{t + e)Bo is within £^ distance Sq of one of the a;"s. By the definition 
(4.9), fx{x'^) >t + e. Any u G Am satisfies 

(4.12) \u - [mx']\ < mdo + d 

for some 1 < i < k. For this same i, 

T(0, [ma;*]) < T(0, u) + T{u, [mx']) . 

T{u, [mx'^]) is stochasticaUy dominated by a sum 5'r'^^ i of [m5i\ i.i.d. Exp(l) ran- 
dom variables, because by (4.12) a direct lattice path from u to [mx'^] takes at most 
mdo + d < m{SQ + d/n) < mSi Exp(l) passage times. Consequently, by (4.10), 
Proposition 4.1, and standard large deviation bounds for exponential random vari- 
ables, for u G Am, 

P{T{0,u) < nt) 
k 

< 5^ P(T(0, [mx']) < m^{x') - msi) + kP{Sl^s^^ > m6i{l + Si)) 

i=l 

< Cikexp{—C2m) 

for some finite constants Ci, C2 > [not the same as those in Proposition 4.1]. By 
substituting this bound in (4.11) we have proved Corollary 4.1 for large enough n. 
By increasing Ci in the statement it follows for all n. D 

Corollary 4.2. For any t, e > there are finite constants Ci = Ci{t,e) > such 
that 

P{n{t - e)Bo ^ B{nt)) < Cx exp(-C2n) 

for all n. 

Proof. The event n{t — e)Bo ^ B{nt) implies that T(0, ti) > nt for some u G 
n(t — e)Bo. The idea for proving this is the same as for the previous Corollary: 
Pick a fine grid of point x^, . . . , x^ G (t — e)Bo. For each u G [n{t — e)Bo\ fl Z^ 
pick the closest [nx^], use the inequality 

T(0, u) < T(0, [nx']) + T{[nx'],u) , 

dominate T([na;*], u) by a sum of i.i.d. Exp(l)'s, note that fx{x^) < t — £, and apply 
Proposition 4.1. We leave the details to the reader. D 

As the last item of this section we prove Proposition 4.1. For this purpose we 
combine a result of Talagrand [30] , stated as the next lemma, with some ideas from 
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Kesten [14] . Suppose {Xi : 1 < z < A^} are independent random variables such that 
< Xi < 1. Let JF be a family of A^-tuples a = (a^ : 1 < i < N) of real numbers. 

Set a = supQ,g;p- ||a||25 where ||a||2 = I Xli=i '^f ] ■ Define the random variable 

TV 

(4.13) Z=supVajX, . 

Lemma 4.2. [30] For any numbers a <h, 

(4.14) P{Z >h)-P{Z <a)< exp [-{h - af/Aa'^] . 

Some comments: We use here the original notation of [30] even though the use 
of Z, a and JF conflicts with our notation elsewhere. But this notation is used only 
for the proof of Proposition 4.1 so it will not cause confusion. Except for Lemma 
4.2 our proof of Proposition 4.1 is self-contained, and the proof of Lemma 4.2 can 
be quickly read from Talagrand's paper [30]: the reader only needs to go through 
section 4.1 which is self-contained, and the proof of Theorem 8.1.1 which relies on 
section 4.1. Our Lemma 4.2 is contained in the proof of Theorem 8.1.1 in [30]. 

Fix X G R'^. To produce the random variable Z of (4.13), we truncate the 
random times of the sites and restrict the set of paths tt. Let X, L > be con- 
stants independent of n, to be chosen below. Let Hi, denote the collection of 
self-avoiding nearest- neighbor paths tt from to \nx\ with length at most Ln: 
TT = {0 = uf* ,w^ , . . . ,w'^~^ .w"^ = [nx]} such that m < Ln. Define truncated 
times by t{u) = t{u) AK for u E Z^. Let T denote the passage time with truncated 
variables and restricted paths: 

T(0, [nx]) = min T(7r) , 

where T(7r) is defined by (4.1) with t replaced by t. 

Let {0 = w^, -u^, w^, . . . ,'U^} be a numbering of the sites that are within £^ 
distance Ln of the origin. Let JF be the subset of { — 1, 0}-valued A^-tuples a = (a^) 
that are negatives of indicator functions of paths tt G Hl'- a G JF iff for some 
n G n^, ctj = — !{«* G tt} for all 1 < z < A^. Define i.i.d. random variables with 
values in [0, 1] by Xi = K~^t{u'^). The random variable Z defined in (4.13) is then 

(4.15) Z = sup f]a.^ = -1 inf Yl ^(^^) = "^^(0' M) • 
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Lemma 4.3. Given e > Q, we can choose the constants L — L{x) and K — K{e, x) 
so that 

(4.16) P{ |f (0, [nx]) - T(0, [nx])\ > ne) < Ci exp(-C2n) 

for all n, for some finite constants Ci, (^2 > 0. 

Proof. T(0, [nx]) has a.s. a unique optimal path vr^ because the site times t{u) have 
a continuous distribution. T(0, [nx]) does not necessarily have a unique optimal 
path, so order the paths in 11^ in some way and let vr^ denote the first path in 
this ordering that is optimal for T(0, [nx]). Divide the estimation in (4.16) into two 
parts: 

P{\f{0,[nx])-T{0,[nx])\ > ne) 

< P{'K^ i Hz.) + P( |f (0, [nx]) - T(0, M)| > ne, 7r„ G H^) . 

The first term of the right-hand side is bounded as follows: 

< P(T(0, [nx]) > n,Lo) + -P( a self-avoiding path tt starts from 

and has length at least Ln but T{t{) < nLo J 
<P(5f[_jl>nLo)+ Yl {2d-l)'-pisl<m-^y 

m^Ln 

Above we used T(0, [nx]) < ^'.^r ^i,, bounded the number of self-avoiding paths from 
with length m by {2d — 1)'^ [at each step at most 2d — 1 sites to choose from], 
and noted that for a path tt of length m, T{ti) is distributed like 5"^, a sum of m 
i.i.d. Exp(l) variables. By choosing Lq = x + 1 and L = L{x) = 2a; + 2 the quantity 
above is bounded by 

oo 

exp(-Cin) + J2 (2c/- l)"'exp(-C2m) < Ciexp(-C2n) 

for constants Ci,C2 > 0. Note that the constants change from one inequality to 
the next. 

To bound the probability P( |T(0, [nx]) — T(0, [nx])\ > ne , 7r„ G 11^) note first 
that if Tin G n^, then necessarily 

< T(0, [nx]) - f (0, [nx]) < J^ [t{w) - t{w)] 
= J2 [tH - K] ■ l{t{w) > K} . 



W&n 
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For paths tv G Hl define the events 

A^ = \Y1 [^(^) ~ ^] ■ 1{^('^) > K}>ne\ and D^ = {tt^ = tt} . 

Let ^TT denote the cr-algebra generated by {t{w) : w ^ tt}. When tt and {t{w) : w ^ 
n} are fixed, A^r is an increasing event and Dyr is a decreasing event of the variables 
{t{w) : w G tt}. We estimate as follows: 

P ( I f (0, [nx] ) - T(0, [nx] ) | > ne , tt^ G n^) 

< P( ^ [t{w) - K] ■ l{t{w) > K}> ne\ 

[ by the FKG inequality ] 

= Y, P{A^)E[P{D^\g^)] 

[ Ajr is independent of ^jr ] 

^ i=l ^ ttETIl 

[ each TT has at most Ln passage times ] 
- ^(X][^(^') - ^] ■ HKu") > K}> ne^ . 

This last probability is < exp(— Cn) if we choose K = K{e, L) large enough so that 

E[{t{u') - K) ■ l{t{u') > K}] < e/{2L). D 

We are ready to prove Proposition 4.1. For convenience, replace e by 8e in the 
statement of the Proposition. Fix e > 0, and choose K, L so that (4.16) holds. 

P{ |T(0, [nx]) - n^{x)\ > 8ne) 
< P ( I T(0, [nx] ) - f (0, [nx] ) | > 4ne) 
+ P{\f{0,[nx]) - nfi{x)\ > 4ne) , 
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SO by Lemma 4.3 it suffices to bound the last probability. By the limit in (4.8) and 
by Lemma 4.3, 

(4.17) P{ |f (0, [nx]) - nij{x)\ > 2ne) < 1/2 

for all large enough n. Recall equality (4.15), use (4.17) and estimate as follows: 

P{ |f (0, [nx]) - n^{x)\ > 4ne) 
= P{ f (0, [nx]) > n^{x) + 4ne) 

+ P(f (0, [nx]) < nfi{x) - 4ne) 
= P{Z < -n^{x)/K - Ane/K) 

+ P{Z>-n^{x)/K + Ane/K) 
<2-P{Z < -nij{x)/K - Ane/K) ■ P{Z > -nij{x)/K - 2ne/K) 

+ 2-P{Z> -nij{x)/K + Ane/K) ■ P{Z < -nij{x)/K + 2ne/K) 

<Aexp[-{2ne/K)'^/ALn] 
= Aexp[-e'^n/{K^L)]. 

We used Lemma 4.2 in the second last step. Each a E J-' has at most Ln entries 
equal to —1 and the rest are zeroes. Consequently a < \J Ln . We have proved 
Proposition 4.1 for large enough n, and it follows for all n by increasing Ci suffi- 
ciently. 

5. The ballistic deposition shape from a seed 

In this section we prove Theorem 1, the almost sure limit for the ballistic depo- 
sition process started from a single seed. Initially 

. X . . r if w = 

(5-1) ^.(0)= .„ ,^ 

Caution to the reader: Z was defined by (4.15) only for the proof of Proposition 4.1. 
Elsewhere in the paper, Z denotes the ballistic deposition process started from a 
seed. The process Z is constructed by the argument of Section 3 on the probability 
space (O, JF, P) of the Poisson jump time processes T = {T^}- Recall that So ^ R'^ 
denotes the closed, convex limiting set for ffist-passage site percolation on Z'^ with 
Exp(l) waiting times. The goal is to prove that for a bounded, positive, concave 
function g defined on the open set intSo, 

(5.2) lim —Z\n^-\{nt) = tg{x/t) a.s. 
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for all X G R'^ and t > such that x/t G intSo- 

The main tool in the proof is the Kesten-Hammersley lemma from subadditive 
ergodic theory. Since Z is unbounded both above and below, we work instead with 
the passage times R{u,h) defined by (4.5). By definition R{u^h) > 0. Since Z 
can always reach cell (w, h) by advancing along each coordinate axis in turn, with 
Exp(l) waiting times, 

(5.3) R{u, h) 4 ^|\.|+^ . 

Here =^ denotes stochastic dominance and 5"^ is a sum of n i.i.d. Exp(l) random 
variables. Recall that \u\ is the l^ norm on Z^. In particular, R{u^ h) has finite 
moments of all order. 

The main property to check is a subadditivity: 

Lemma 5.1. For u,v ^TJ^ and h,k E Z_|_, there exists a random variable R{v, k) 
such that 

(5.4) R{u + v,h + k) <R{u,h) + R{v,k), 

and R{v, k) is independent of R{u, h), and equal in distribution to R{v, k). 

Proof. To define R{v, k), we start a new Z-process at time R{u, h) from a seed in 
cell {u,h). This new process Z is dominated by the original Z-process, hence the 
inequality (5.4). 

For the reader not familiar with these types of arguments, here are the details: 
Think of T = {T^ : ty G Z*^} as an infinite- dimensional vector of Poisson point 
processes T"", indexed by time t G (0, oo). Let J-'t be the u-algebra generated by 
the restriction of T to the time-interval (0,t]. Then R{u, h) is a stopping time for 
the filtration (J-'t)- Restart the Poisson processes at time R{u, h) and translate the 
index by u to get new point processes T = {T'^ : w G Z'^} where T^ = |^7"^+"' — 
R{u, h)~\ n (0, oo). (The subtraction means that epochs of 7"""'""' are translated back 
R{u, h) time units.) 

Sublemma. T is a collection of i.i.d. rate 1 Poisson processes on (0, cxd), inde- 
pendent of R{u,h). 

Proof of sublemma. We ignore the spatial translation w ^^ u + w. R = R{u, h) is 
the stopping time, and Tr is the cr-algebra up to time R, consisting of events A such 
that An{R<t} e J^t for aU t. Topologize the space of families T = {T™ : ty G Z^} 
of point processes on (0, cxd) by the product of the vague topology on point processes. 
Let (7 be a bounded continuous function on that space. The claim is that for any 
A G JFjj and any such (7, 

(5.5) E[lAg{f)]=P{A)E[g{T)]. 
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Since R itself is jFjj-measurable, the lemma follows. 

Take a sequence of discrete stopping times Rn that decrease down to R, almost 
surely. For example, Rn = E/t k2-''l{{k - 1)2"'^ < R< /c2-'^}. Let f^ denote T 
translated back R^ time units and restricted to (0, oo). A G J-'r implies A G J-'r„, 
and we get: 

oo 

E[1a 9{fn)] = J2 E[^A l{Rn = U}g{fn)] 

i=l 

CXD 

= J2E[lAHRn = U}g{{T - U) n (0, oo))] 



1=1 

CXD 



J2E[lAHRn = U}]E[g{r)] 



1=1 
= P{A)E[g{T)] . 

The {ti} are the possible values of Rn- The point of the above calculation is of 
course that, when ti is a deterministic number, A fl {Rn = ti} G J-'t^, while T — ti 
restricted to (0, oo) is an i.i.d. collection of Poisson point processes independent of 
J-'t^. Now let n -^ oo, and observe that Tn ^> T a.s. so that £'[1^(7(7^)] converges 
to i?[lAfi'(^)], and we obtain (5.5) in the limit. D 

Return to the proof of Lemma 5.1. Define the process Z as a function of T, 
exactly as Z is a function of T, with an initial seed at the origin: Z^{Q) = —00 ■ 
l{w 7^ 0} as in (5.1) for Z. Let R{v, k) be the time when Z gets at or above cell 
{Vj k). Then R{vj k) is independent oiR{u, h), and has exactly the same distribution 
as R{v,k) defined by (4.5). 

Consider the processes a and p defined by 

ayj{t) = Zyj{R{u, h) + t) and pw{t) = h + Zu,-u{t). 

Then pwiS^) = —00 except at w = w, where 

p„(0) = h + Zo(0) = h = Z^{R{u, h)) = (7,(0). 

The jump times of a^^ are given by T^ — R{u^h), and those of pu, by T^~^ = 
T"" — R{u, h). So processes a and p use the same Poisson jump times and initially 
f7(0) > p(0). By Lemma 3.2 a{t) > p{t) for all t. Take t = R{v, k). Then 

Zu+v{R{u, h) + R{v, k)) = au+v{R{v, k)) 

> pu+v{R{v, k)) = h + Zy{R{v, k)) 
= h + k, 
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which says that Z^^^ has reached height /i + /c by time i?(w, K) + R{y, k). This 
imphes (5.4) and completes the proof of Lemma 5.1. D 

Inequahty (5.4), the existence of moments, and the Kesten-Hammersley lemma 
as given on p. 20 of [28] imply that for all u E Z'^ and h G Z_(_ there exists a number 
7(tt, h) such that, for any positive integer m, 

(5.6) lim —. — R{2^mu,2^mh) = '^{Ujh) a.s. 

The first task is to extend (5.6) to a genuine limit, and for that we use the 
following continuity property of the passage times. It is an immediate consequence 
of (5.3), (5.4) and the monotonicity of R{u, h) in the /i- variable. 

Lemma 5.2. For any u,v E Z*^ and h,k E Z_|_, 

(5.7) Riu,h)-Riv,k)4Sl^_,^^^^_,^+. 

Now regard (tt, h) fixed, and also fix e > and an integer m large enough 
relative to e and (w, h) [how large m needs to be is seen shortly] . For large enough 
n, pick first j = j{n) so that 2^m < n < 2^~^^m. Then it is possible to find a 
k = k{n) G {0, 1, . . . , m — 1} such that 

(5.8) 2^{m + k)<n<2^{m + k + l). 

We write 

2^m + k) 1 



(5.9) 



n 2^{m + k) 



R{2^ {m + k)u,2^ {m + k)h) 



< —R{nu^nh) H — U{n) 
n n 



where the error 

U{n) = {R{2^ {m + k)u,2^ {m + k)h) - R{nu,nh)} y 4 S^j^ 

for a constant L — L{u,h), by (5.7). The Cramer rate function for Exp(l) is 
k{x) = X — 1 — logx, so we get the estimate 

P{U{n)>ne)<p(^Sl,^>2^L^^ 

(5.10) <p(sl,L>2^Lme/L\ 

<ex.p[-2^LK{me/L)] 
< exp[— Cn] , 
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where C > is a constant. In the above calculation we used 2^m < n < 2^^^m 
and k{x) > x/2 for large enough x, and took m large enough. 

Let n ^ oo in (5.9), so that j ^ oo also. Even though k = k{n) varies with n, 
it has only finitely many possible values so the limit (5.6) happens on the left-hand 
side of (5.9). Note that 

2^{m-\-k) 171 + k m 



n m + k + 1 m + 1 

The error n~^U{n) vanishes a.s. by the estimate (5.10) and Borel-Cantelli. We get 

1 TTl 

lim inf —i?(n'U, n/i) > 7(w, /i) a.s. 

n^oo n m + 1 

A similar argument works for the linisup. Let m ^ oo, and we have this inter- 
mediate statement: for all w G Z'^ and h G Z+ there exists a number 7(^5 h) such 
that 

(5.11) lim —R{nu^nh)='-y{u,h) a.s. 

n— >oo 77, 

From (5.7) we get a Lipschitz property for 7: 

(5.12) |7('u, h) - 7(1), k)\<\u-v\ + \h-k\, 
while (5.11) gives homogeneity: for nonnegative integers m 

(5.13) '-yimu.mh) = m'-y{u,h) . 

(5.13) permits us to define 7 unambiguously for (x,6) G Q^ x Q_|_ by 

1 

(5.14) 7(x, 6) = — '-y(mx,mb) 

m 

where m is any positive integer such that (mx^mb) G Z^ x Z_|_. By an estimation 
similar to that used in (5.9), the limit in (5.11) can be extended to all (w, /i) G 
Q^ X Q+. 

The final step is to extend the limit in (5.11) so that it holds outside a single 
exceptional P-nuU set for all (-u, h) = (xo,6o) ^ R-"^ x R'+- The Lipschitz property 
(5.12) continues to hold on Q'^ x Q_|_ for the extension of 7 defined by (5.13), so 
we can extend 7 uniquely to a Lipschitz function on R'^ x R^. For rational (x, 6), 
integers n, and rational S > define the event An^5{x, b) by 

An^sixj b) = < there exists {v, /c) G Z'^ x Z+ such that 

(5.15) I [nx] — v\ + \ [nb] — k\ < n5 but 

|i?([na;], [nb]) - R{v, k)\ > 2nd 
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By (5.7) and by standard large deviation bounds for exponential random variables, 

P(A„,5(x,6)) < Cin^+iexp[-C2n] 

for finite constants Ci = Ci{x, b, 6) > 0. This bound is suitable for Borel-Cantelli. 
Thus at this stage the following holds with probability 1: for each (x, 6) G Q'^ x Q_|_ 
and rational 5 > 0, (5.11) holds with (w, h) = (x, 6), and for large enough n, 

(5.16) \R{[nx], [nb]) - R{v, k)\ < 2nd 
for all cells (f , k) that satisfy 

(5.17) I [nx] — v\-\-\ [nb] — k\ < nS. 

Now let (xojbo) G R'^ x R_|_. Choose rational (x, 6) such that |x — xqI + |6 — 6o| < 
d/2. Take {v,k) = {[nxo],[nbo]). Then (5.17) holds for large enough n, and by 
letting n — *> oo in (5.16) we get 

7(x, b) — 26 < liminf —R{[nxo], [nbo]) < linisup —R{[nxo], [nbo]) < 7(x, b) + 26. 

Take rational (x, b) that converge to (xq, bo), and use the Lipschitz continuity of 7. 
This proves the limit in this proposition: 

Proposition 5.1. There is a homogeneous, subadditive, convex Lipschitz function 
7 on H'^ X R_|_ such that, outside an event of probability zero, 

(5.18) lim -R{[nx], [nb]) = 7(0:, b) 

n^oo n 

for all {x,b) G R"^ x R+. 

To prove Proposition 5.1, it remains to argue the properties of 7: The subad- 
ditivity (5.4) implies the corresponding subadditivity for 7. This subadditivity is 
preserved by the extensions of 7 to rational and real points of R'^ x R+. Same 
is true of the homogeneity (5.13). Homogeneity and subadditivity together imply 
convexity. 

At this point we have not ruled out the possibility that 7 = 0. (5.3) gives the 
upper bound 

(5.19) -f{x,b) <\x\ + b for all (x,6) G R"' X R+. 

To get a positive lower bound for 7, we construct another growth process in Z'^ x Z_|_ 
whose height dominates Z, and that has a simple structure so that its spread is 
easier to bound. 
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Instead of just focusing on the height Zu{t) of the growing cluster, let us denote 
by Z{t) the actual set of occupied cells at time t. Z{t) is a subset of Z*^ x Z_|_, and 
initially Z{Q) = {(0,0)}. The rule of evolution for Z{-) is this: At epochs of T", 
the top growth cell above u is annexed to Z. By definition, the top growth cell 
above site u is (w, /c) G Z'^ x Z_(_ with maximal k subject to the condition that 

(5.20) (w, k) ^ Z, but either {u^k — 1) E Z or {u + z,k) E Z for some z G N . 

If no finite k satisfies this condition, nothing is annexed to Z. 

Define a different growing cluster W by stipulating that all growth cells (and not 
just the top one) join independently with rate 1. A cell (-u, k) is a growth cell for W 
if (5.20) holds with W instead of Z. So note in particular that the cluster W grows 
"sideways" and "up" in Z<^ x Z+, but not down. Set initially >V(0) = {(0,0)}. 

To construct W we employ a collection {T^ : u G Z'^, q G Z_|_} of i.i.d. Poisson 
point processes. The nonnegative integer q labels the growth cells from top down, 
so that the top growth cell above site u is assigned label q = 0, the next highest 
growth cell gets label q = 1, and so on. More formally, we can define 

gu,o = max{/c : (w, k) is a growth cell for W} 

("(/" for growth) and inductively for q>l 

gu,q = max{0 < k < gu,q-i '■ ("W, k) is a growth cell for W} . 

By convention, the maximum of an empty set is — cxd, so if the current cluster W 
has exactly m growth cells above site u, then gu,q = — oo for q > m. The gu,q^s are 
of course functions of time too. The precise rule of evolution for W is this: 

if r is an epoch of T^ and gu,q > — oo, then >V(r) = W{r—) U {{u,gu,q)}- 
Update the (7„,g's forv = u and v E u + M. 

The argument given for the construction of ballistic deposition in Section 3 
does not work for W because each site u now has infinitely many Poisson processes 
{Tq '■ Q ^ Z_|_} attached to it. However, we can easily see that, given any finite time 
ti, at most finitely many TT'-processes are involved in constructing the dynamics up 
to time ti: Starting with the seed at (0, 0), let ti, T2, T3, . . . be the successive waiting 
times for adding the second, third, fourth,... cell to the existing cluster. Since each 
new particle adds at most 2d + 1 growth cells, Tn is stochastically larger than an 
Exp(2(in + n) random variable. Consequently '^^r^ = 00 a.s. and only finitely 
many steps are needed (and only finitely many Poisson processes T^ inspected) to 
construct W{t) for < t < ti. 
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Now couple W(t) and Z{t) by letting {Tf^} be the Poisson processes that govern 
the evolution of Z. In other words, both Z and W annex their top growth cell above 
u simultaneously at epochs of 7^". At epochs of T^ for g > 1 the Z-process does 
nothing, while W may add other growth cells as stipulated in rule (5.21). Since the 
top growth cells are not necessarily the same for Z and W it does not follow that 
W(t) contains Z{t), but it does follow that the height of W always dominates the 
height of Z. If we let 

Wuit) = max{A; >0:{u,k)e W{t)} 

when some (w, k) lies in W(t) and Wu{t) = — oo otherwise, we get this inequality: 

(5.22) Zu{t) < Wu{t) for all u e Z'^ and t > 0, a.s. 

Proof of (5.22) is by induction on jumps [which are only finitely many in any 
finite time interval, almost surely, as argued above]. As long as (5.22) holds, no 
top growth cell of Z can be above the corresponding top growth cell of W, and 
consequently the next jump preserves (5.22). 

To make use of (5.22) we redefine W as a first-passage problem. Give the cells 
i.i.d. Exp(l) random waiting times {t{u,h) : {u,h) G Z^ x Z_|_}. Consider self- 
avoiding nearest-neighbor lattice paths tt = {(i'°, k^), {v^, k^), • • • , (f"^, k^)} whose 
admissible steps are of these types: for each i = 1, . . . , m, 

(5.23) {v\ k') - {v'-^,k'-^) = (0, 1) or (icp, 0) for some p = 1, . . . , rf, 

where ei, . . . , e^ are the d standard basis vectors in R*^. In other words, inside a 
layer Z*^ x {/c} an admissible path it takes arbitrary nearest- neighbor steps subject 
to self-avoidance, and across the layers n moves only up, not down. The admissible 
steps are chosen to match (5.20). The passage time of such a path tt is M(n) = 
Xli^i ^(^^5 ^*)- The passage time of cell (ti, h) is 

M{u,h) = infM(7r) 

TT 

where the infimum ranges over paths n of the above type from (0, 0) = {v^, k^) to 
(v^,^™) = {u,h). Let 

W{t) = {{u,h) :M{u,h) <t} 

be the growing cluster. For W we get a bound easily by counting self-avoiding 
paths. 
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Lemma 5.3. For any finite constant j3, there exists a positive v = z/(/3) such that 
for any finite integer K , 

(5.24) P { W(t) contains a cell {u, h) such that \u\+h>K\ < 2 exp(-ii:/3) 

as long as t < Ku. 

Proof. Let w denote an admissible path fixed to start at the origin (0, 0) = (f^, /c°). 
The probabihty in (5.24) is at most 

PI there exists a tt through at least K cells with M(7r) < t 



oo 

< 



y P[ there exists a tt through exactly j cells with M(7r) < t 

j=K ^ 



< J2{2d + iyp{s]<t) 

j=K 

< Y,expl-j[K{t/K) -log{2d+l)]\ 



j=K 

oo 

< J]exp(-/3j)<2exp(-K/3). 

j=K 

We used again the Cramer rate function k{x) = x — 1 — logx for the Exp(l) 
distribution. The inequality PiSj < t) < exp[—JK,{t/j)] is valid for t < j. Take 
t/K < V with V small enough so that k{v) > /3 + log(2(i + 1) + 1. D 

Lemma 5.4. The processes W(-) and W(-) are equal in distribution. 

Proof. We can regard both processes as jump processes on the countable state 
space of finite connected subsets of Z'^ x Z+ that contain the origin. Both start 
from {(0,0)}. Both processes add new admissible cells independently at rate 1. 
Comparison of (5.20) and (5.23) shows that an admissible new cell, or growth 
cell, is the same for both W(-) and W(-). Thus the two processes have identical 
infinitesimal rates. D 

Combining (5.22) and Lemmas 5.3 and 5.4, we get a lower bound for 7: 

Lemma 5.5. There exists a positive constant v such that 7(x,6) > i'(\x\ + 6) for 
all {x,b) G R^ X R+. 
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Proof. Pick u < u{l) = the constant given by Lemma 5.3 for (3 = 1, and set 

t = z/(|a;| + b). Then by Lemma 5.3 for K = \[nx]\ + [nb], 



P( R{[nxl [nb]) <nt) < P( Z[„,](nt) > [nb] 

<plW[^.]{nt)>[nb] 
< 2exp(— |[nx]| — [nb]). 



Thus R([nx], [nb]) > nt for large enough n, a.s. D 

Now we have 7 bounded both above and below. Finally, we convert the limit in 
Proposition 5.1 to that of Theorem 1. We need one more property for 7: 

Lemma 5.6. For x G intSo; there is a unique finite h > such that 7(0:, h) = 1. 

Proof. For this proof we connect ballistic deposition on Z*^ x Z+ with first-passage 
site percolation on Z*^. Recall the definition of the limit n{x) in (4.8). By (4.7) 
and Proposition 5.1, n{x) = 7(x,0). By definition, Bq = {x : n{x) < 1}, so by 
homogeneity ^{x) < 1 for x G int Bq. Hence for such x also 7(0;, 0) < 1. By Lemma 
5.5 and the continuity of 7 there is some h > such that 7(0:, h) = 1, so it remains 
to rule out the possibility of having < ho < hi such that 7(x, ho) = 7(x, hi) = 
1. But this and convexity would imply 7(^,0) > 1, contradicting what was just 
concluded. D 

Proof of Theorem 1. By the previous lemma, a positive function g on intSo is 
uniquely defined by the equation '^{x,g{x)) = 1. The lower bound of Lemma 5.5 
gives an upper bound for g, and g is concave by the convexity of 7. A finite 
concave function on an open convex set is continuous by Theorem 10.1 in [Rf]. By 
the homogeneity of 7, 

(5.25) 'y{x,tg{x/t)) =t. 

By the uniqueness in Lemma 5.6 and the monotonicity of 7, /i > tg{x/t) [h < 
tg{x/t)] implies 7(x, h) > t [7(0:, h) < t]. 

Fix X G R'^ and t > so that x/t G int So- Let So > 0, and pick Si G (0, So). Set 
b = tg{x/t) — Si, and pick e > so that 7(0:, b) < t — e. Fix a large number mo so 
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that mo(5o — 5i) > 1 to take care of the effects of integer parts. 



hminf n "'^Z'r^^i(nt) < tg{x/t) — 5o 

OO OO f s 

^ n U { %-] (^^) < [^i9{x/t) - n5i] I 
C Pi y J R{[nx], [nb]) > w^ix, h) +ne\ 

m=mo n=m ^ ^ 

C I hm sup n~"'^i?([na;], [nb]) > 7(x, b) + e> . 

l_ n — >oo J 

The last event has probabihty zero by Proposition 5.1. Similarly we show that 

\iBisupn~^ Z^rixjin-t) < tg{x/t) a.s., 

n — »oo 

and the limit in Theorem 1 is proved. The invariances of g follow from the corre- 
sponding invar iances in the distribution of Z. D 

6. The shape from a translated seed 

This section is a purely technical extension of the limit from a seed (Theorem 1) 
derived in the previous section. The proof of the hydrodynamic limit (Theorem 2) 
uses the coupling (3.9) which forces us to consider simultaneously the whole family 
{Z'"} of processes. Recall that Z^ stands for the process that grows from a seed in 
cell {v, 0), as defined by (3.6). We need a limit where the initial seed and the initial 
time point are translated as the limit is taken. 

In Section 3 we defined the processes {Z^} on the probability space (O, JF, P) of 
the Poisson jump time processes T = {T^}. Now we extend the Poisson processes 
to the entire real line (—00,00). Since Poisson points on (— oo,0] and (0, 00) are 
independent, this is the same as replacing the original probability space (O, JF, P) 
with a product space 

{n,T,p) = (0° X o, j^o ® j^, p° ® p) 

where (0'^,JF°,P°) is the probability space of the Poisson processes on (— oo,0]. 
A sample point of the product space is cJ = (a;'^,a;), where u still represents a 
realization of the i.i.d. collection T = {T" : u G Z'^} of Poisson point processes 
on (0,00), while lj^ represents a realization of these processes on (— oo,0]. The 
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processes Z"" are defined on the space Q in the obvious way, by ignoring the uP- 
component: Z^(t)(c^°, w) = Zl{t){u)). _ 

Let Qs for — oo < s < oo denote a time translation on the space O, so that the 
epochs of d^ are those of U translated s time units backward. (In other words, 
around the time origin d^^ looks like To around time point s.) The random variable 
Z^{t) o 9s is computed by reading Poisson jump times from time point s onwards 
and by letting the ballistic deposition evolve for a duration t. The passage times 
are defined as before: 

R^u, h)oes= inf{t >Q:Zl{t)oes>h}. 

Notice that R'"{u,h) o 0^ is the amount of time it takes to grow up to height /i, 
and not the time point of the Poisson processes when this happens. Note further 
that if s < then the random variable Z^ (t) o Og uses the Poisson processes also 
for negative times. 

Theorem 4. Let g be the limiting function of Theorem 1. Then the following holds 
with probability 1: For all x,y E R'^, s G R, and t > such that y E x + t ■ int Bq, 

(6.1) hm -Z\:%t)oe^s=tg(''~y 



n^oo n '■^^' \ t 

The remainder of this section proves Theorem 4. We construct step by step an 
event F on O that satisfies -P(r) = 1 and on which the convergence (6.1) holds. 
As in Section 5, we prove the theorem by proving the convergence of the passage 
times: 

(6.2) lim -Rl'^y^ ([nx], [nb]) o O^s = l{x - y, b) 

n^oo n 

for all y G R'^, (x, b) G R'^ x R_|_, and real s. 

The reader who is willing to accept the result may ignore the proof without 
loss of continuity and proceed to the next section where the hydrodynamic limit is 
proved. 

Step 1. First we define F as the event on which 

(6.3) lim -^i?2^"^« (2^mw^ 2^m/i) o O^^ms, = l{u - v, h) 

for all V G Z'^, [u, h) E Z'^ x Z_|_, rational Sq, and all m G N. 

We need to argue that -P(F) = 1. The proof of the Kesten-Hammersley theo- 
rem on p. 20-23 of [28] shows that convergence along powers of 2 only depends 
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on the distributions of the random variables. Let F^ denote the distribution of 
R^^{nu,nh) o 6ns^y Then Fn is also the distribution of R{n{u — v),nh), and in- 
equality (5.4) shows Fm+n > Fm-^Fn- This, and the existence of second moments, 
is what is needed for the almost sure convergence in (6.3). 

Step 2. Now define F to be the event where the requirement of Step 1 holds, 
and in addition properties (6.4)-(6.5) below, which are to hold for all y G Q'^, 
(x, h) G Q'^ X Q_|_, rational Sq, and rational Si, 5 > 0: 

(6.4) For large enough n, \B}-'^y^ {[nx], [n6]) — i^t"-?/] (^^ /j) | < 2n5 for all cells {v, k) 
that satisfy | [nx] — v\ + \ [nh] — k\ < n5. 

(6.5) For si > there exists 5o(si) > such that if 5 < 5o(si), then for large 
enough n, Zv'{nsi) o 6nso > for all sites v such that 1^; — [ny]\ < nS. 

Still -P(F) = 1. We already argued around (5.15)-(5.17) that (6.4) can be satis- 
fied almost surely. By Lemma 4.1, condition (6.5) can be viewed as a percolation 
question: For each n define a first-passage percolation cluster centered at [ny] by 

B{t) = {vez'':z^;'y\t)oenso>o}. 

If 5o('Si) is chosen so that \v — [ny]\ < ndo{si) implies v G [ny] + n(si/2)i?o, then 

(6.5) holds a.s. by Corollary 4.2 and the Borel-Cantelli lemma. 

The next step is to improve the convergence in (6.3) to a genuine limit on the 
event F. Fix m for the moment. As in (5.8), for large enough n there are j = j{n) 
and k = k{n) G {0, 1, . . . , to — 1} such that 

(6.6) q{n) = 2^(m + k) < n < 2^m + k + l). 

Keeping v G Z'^, (tt, /i) G Z'^ x Z+ and so ^ Q fixed, we use properties (6.4) and 
(6.5) of F to write 

i?^('^>(Q(n)w,Q(n)/i)o^^(„)(,^_,) 

(6.7) < (n - q{n))so + q{n)s + i?"" {q{n)u, q{n)h) o 6'„,^ 
< (n — q{n))sQ +ns + R^^ {nu, nh) o 6*^50 + 25n . 

This argument will be used several times, so we go over it once carefully. The first 
inequality in (6.7) is valid on the event F for a certain s = s{ni) > that satisfies 
lim^^oo s{m) = 0, and for all large enough n, by property (6.5) for this reason: The 
ballistic deposition process Z'^^'^'^ o 6'g(^)(sg_s) starts from a seed in cell (g(n)t', 0) 
when the Poisson process clock is at q{n){sQ — s). One way for this process to reach 
cell {q{n)u, q{n)K) is to first spend at most time (n — q{n))sQ + q{n)s to reach cell 
(nf , 0), and from there follow a new process Z'^'" oOns^ that starts when the Poisson 
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process clock is at usq. Since < n — q{n) < n/m, there is a 5 = 6{m) > such 
that hm^^oo S{m) = and \q{n)v — nv\ < nd if n is large enough. Hence by (6.5) 
we can choose s = s{m) > such that Z'^^'^')'" o 6'g('„)(s_SQ) reaches cell {nv,0) in 
time q{n)s if n is large enough. 

The second inequality in (6.7) follows from (6.4), again for a, 6 = d{m) > such 
that lim^^oo S{m) = 0, and for large enough n. Also q{n) < n was used. A note 
about terminology: We say that Z reaches cell (w, h) in time t if Zu{t) > h. 

Let n ^ oo in (6.7), use the limit (6.3) and that n — q{n) < n/m to conclude 
that on r, 

(6.8) liminf -i?""(nw, nh) o 6'„,„ > -^{u - v,h) - is + 25) - — . 

n^oo n m m 

Now let m -^ oo, and concurrently we can take s — > and 6^0. 
To handle the limsup we use the same argument in a reverse way: 

R''^{nu,nh)oer,so 

(6.9) < iq{n) -n)so + q{n)s + W^'^^ [nu, nh) o ^,(n)(.+so) 
<ns + i?^(")^ {q{n)u, q{n)h) o ^,(^)(,+,„) + 25q{n) . 

Now the thinking goes like this: The process Z^'" o 0^^^ starts from a seed in cell 
(nv, 0) when the Poisson process clock is at usq. One way for this process to 
reach cell {nu, nh) is to first spend at most time (g(n) — n)so + q{n)s to reach cell 
(g(n)i', 0), and from there follow a new process Z'^^'^^'" o 6'q(^)(s_|_sy) that starts when 
the Poisson process clock is at q{n){s + sq). Again, there is a 5 = S{ni) > such 
that lim^^oo 5("^) = and \q{n)v — nv\ < nd if n is large enough. By (6.5) we 
can choose s = s{ni) > such that lim^n^oo s(?n) = and Z^'" o 6*^50 reaches cell 
{q{n)v, 0) in time {q{n) — n)so + q{n)s if n is large enough. The second inequality 
in (6.9) follows from (6.4). 

Let n ^ 00 in (6.9) to conclude that on the event F, 

(6.10) limsnp -R'''' {nu,nh) o 9nso < -f{u - v,h) + {s + 26) . 

Let TO ^ 00, and concurrently we can take s — * and 5^0. 

The limits (6.8) and (6.10) permit us to strengthen the definition of F, without 
losing P(F) = 1: 

Step 3. The requirements of Step 2 hold on F, and also 

(6.11) lim -R^"" inu, nh) o O^s,, = 7(w - v, h) 

n— >oo n 
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for all V e Z^, {u, /i) G Z^ X Z+, and sq G Q. 

As previously in Section 5, next we extend the limit in (6.11) to all rational sites 
and cells. Fix y E Q'^, (x, b) E Q'^ x Q_|_, and Sq E Q. Fix an integer m such that 
my E Z'^ and {mx.mh) G Z'^ x Z_|_. For each n, pick j = j{n) so that 



Then on F, 



lim 



q(n) = jm < n < {j + l)m. 



^m(-)y{q{n)x,q{n)b)o9^(^^^ 



n^oo q[n) / IV ; u 



1 

(6.12) j^oo jm 



lim ^R>'^y{jmx,imh) oOjmso 



= — '-yimx — my, mb) 
m 

= l{x-y,h) 

by the homogeneity of 7. The inequalities (6.7) and (6.9) can now be repeated, by 
replacing v,u,h by y^x^b, and by inserting integer parts where appropriate: [ny], 
[nx], [nb]. 

We conclude that -P(F) = 1 for the event F defined as follows: 

Step 4- The requirements of Step 2 hold on F, and also 
(6.13) lim -i?["^l (M, [nb]) o e^so = li^ - y, b) 

for all y E Q'^, (x^b) E Q'^ x Q_|_, and rational Sq. 

The last extension is to arbitrary y E R'^, (x, 6) G R'^ x R_|_, and real s. Let 
£ > be rational, and pick rational d E (0, e/2) so that also d < 5o(£/2) [5o(£/2) as 
defined in (6.5)]. Pick rational yo, xq, bo, sq so that 

\yo -y\ + \xo - x\ + \bo -b\ + \so - s\ < 5/2 . 

Once more we repeat the reasoning that justified (6.7) and (6.9) to write 

<i?["^o](M,[n6])o^^,„_,)+25n 

< n{s -so + e)+ i?[^^l {[nx], [nb]) o O^s + 26n , 
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and 

< n{so -s + e)+ i?["^"] (M, [nb]) o ^„(,,,+,) 

< n(so - s + £) + i?'"^"' {[nxo], [nbo]) o ^n(so+e) + 25n . 
Let n — !► oo to get 

7(a;o - yo, ^o) - 25 - (s - sq + e) 

< limmf-i?["^](M, [n6]) o6'ns < limsup -i?['^^](M, \nb])oens 

< 7(0:0 - yo, ^0) +26 + {so- s + e). 

As we let yo, 3:0, ^0, sq approach y, x, 6, s we can take 5, £ ^ 0. This proves that the 
hmit (6.2) holds on the event F, and concludes to proof of Theorem 4. 

7. The hydrodynamic limit 

In this section we prove Theorem 2. First we work under Assumption A. The 
processes a^ are constructed according to the description of Section 3 on a probabil- 
ity space (n, JF, P) on which are defined the Poisson jump time processes T = {T"} 
and, statistically independently of T, the sequence of initial interfaces {cr'^(O)}. All 
processes a" use one and the same realization of the Poisson jump time processes 
{T"}. By Corollary 3.1 there is a single version of the family of ballistic deposition 
processes {Z""} grown from seeds as defined by (3.6) that satisfy 

(7.1) <(t)= sup{<(0) + Z:(t)} 

for all n G N, tt G Z^ and t > 0, almost surely. Recall that inside the braces the 
correct convention is 00 + (—00) = — 00. The goal is to prove that on some event 
of full probability, for all x G R'^ and t > 0, 

(7.2) hm la^ {nt)=^{x,t). 

For technical reasons we extend the Poisson processes to the entire real line 
(—00, 00) as was done in Section 6. Let (0°, JF*^, P°) be the probability space of the 
Poisson processes on (— oo,0], and consider the product space 

(0,7, P) =(0° xO,J^o®J^,pO®P). 
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A sample point of the product space O is U = (uj^^uj), where w G O represents a 
reahzation of (T, {o""^(0)}) as above, while uj^ G O*^ represents a realization of the 
Poisson processes on the negative time line (—00, 0]. 

If we can prove (7.2) in the product space O, it follows for the original space O 
too. For suppose F C O is an event such that -P(F) = 1 and (7.2) holds on F. Pick 
u)^ G Q'^ such that the a;*^-section 

r^o = {c^gO: (uj^.uj) gF} 

satisfies P(F^o) = 1. Then (7.2) holds on the full-probability event F^^o because 

(7.2) depends on uj only, and not on uj^. So for the remainder of this section assume 
that all Poisson jump time processes are defined for all real times, and we are on 
the event of full probability where Assumption A and Theorem 4 from Section 6 
are valid. 

Recall that ip is defined for (x, t) G R'^ x (0, 00) by 

(7.3) t(;{x,t) = sup {iJo{y)+tg{{x-y)/t)}, 

y6x+t-int Bq 

and on R'^ x {0} set ^/'(^^jO) = V'o(3^)- This first lemma is a consequence of the 
continuity and boundedness of g on intSo- 

Lemma 7.1. Assume i/jq is a continuous [—00, +oo]-valued function on H'^. Then 
ifj is a continuous [—00, +oo]-valued function on R'^ x [0, 00). 

Proof. Let {xmtn) -^ {x,t) in R'^ x (0, 00). We first argue 

(7.4) Yim.iniip{xn,tn)>ip{x,t). 

n — >oo 

Suppose i/j{x,t) > —00 [otherwise (7.4) holds trivially]. Let M < 'i/j{xjt), and pick 
y E X + t ■ int Bq so that 

My)+t9{{x-y)/t) >M. 

Since y E Xn +tn ■ int Bq for large enough n, 

liminf ilj{xn,tn) > liminf<^ i/joiy) +tng( -^ 

n — »oo n — >cxD I \ J^ 

= My)+tg{{x-y)/t) 

>M. 
The equality above follows from the continuity of g on intSo- This proves (7.4). 
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To show 

(7.5) Mo = limsupV'(a::n,^n) < ^(a:^,^) , 

n — >oo 

pick a subsequence Uj so that liuij^oo i^i^nj ^tn.) = Mq. We may assume Mq > — oo 
and il:{x,t) < oo. Let M < Mq, and for large enough j find yj G Xuj + tn. • intSo 
so that 

'0o(2/j) +t^^g{{x^^ - yj)/t^j) > M . 

Since Bq is compact and {xn.,tn.) -^ (ic,t), we may pass to a further subsequence 
so that yj — > y. Define 

y'j =x-—{xn, -yj). 

Then y'j — * y, {xuj — yj)/tnj = {x — y'j)/t, and y'^ G x + 1 ■ int Bq. We wish to argue 
that, for large enough j, 

(7.6) -oo < %{yj) , -iAod/j) < oo . 

By continuity of V'Oj this will follow from showing that — oo < i(^o{y) < oo. We have 



and 



i^oiy) = hm iJoiy'j) < i^{x,t) < oo 



i/joiy) = lim i^oiyj) > M - t\\g\\^ > -oo . 



(7.6) is verified. Now for j large enough, 

M < ipoiyj) +tn^g{{xn, - yj)/tn,) 

= My'j) + tgiix - yj)/t) + {t^^ - t)g{{x - y'^)/t) 

+ i^o{yj) -i^oiy'j) 
< ijj{x, t) + |tn, - 1\ ■ \\g\\^ + ^o(?/j) - i^o{y'j) ■ 

(7.6) was needed to have the difference i^oiyj) — i^oiy'j) defined and convergent to 
0. Letting j /" oo and M /" Mq gives (7.5). This proves continuity on R'^ x (0, oo). 
We omit the similar but shorter argument for the case (a;^, t^) — > (x, 0). D 

One half of the goal (7.2) follows immediately. For arbitrary y G x + t ■ intSo 
such that y G Iq, set u = [nx], v = [ny] and replace t by nt in (7.1), and use 
assumption (2.5) and Theorem 4 from Section 6 to get 

^^•^^ > liminf|laf;^](0) + -^[5{(nt)| 

^My)+t9{{x-y)/t). 
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Note that even though the random variable inside the braces may equal oo + (— oo) 
for finitely many n, eventually n~^Zl^^Unt) is finite because its limit tg[{x — y)/t) 
is finite. For each fixed (x,t) take supremum over these admissible y's to get, by 
the continuity of the functions involved and by the denseness of Yq, 

(7.8) lini^infia[;,](nt)>^(x,t). 

(7.8) holds simultaneously for all (x,t) outside a single event of zero probability. 
The converse 

(7.9) limsup-(7f,,](nt)<^(x,t) 

is where the work is. First we reduce the problem to rational {x,t). Suppose (7.9) 
holds almost surely for all (x, t) E Q'^ x [Qfl (0, oo)] . For rational (x, t) and rational 
s > 0, consider the event 

Dn^s{x, t) = [a^{nt) < o"?^ i(?it + ns) for all sites v 

such that 1^; — [nx]\ < ns/2^ . 

To bound the probability of the complement D'^ s(^'^)' ^^ ^5 {^1^)1 ■s? cind v such 
that \v — [nx] \ < ns/2. Set k — a^^^nt). For times r > nt, define a new process p by 
p(r) = k + Z^{r — nt)o6nt- By the monotonicity lemma 3.2 /9(nt+ns) < a'^{nt+ns)^ 
so cr^t^-iint + ns) < k implies that P[nx](^^ + ^■s) < k which in turn is equivalent to 
Z|^^i(?is) o 9^f = — cx). We can estimate as follows: 

P(D:,,(a:,t)) < Yl Pi<(^t) > cr]l,4nt + ns)) 

v:\v — [nx] |<ns/2 

< E P{Zln4^s) = -00) 
v:\v — [nx] |<ns/2 

< Cin^P(5[L/2] > ns) 
<Cin^exp(-C2?i). 

We used the fact that the time for process Z"" to cover cell ( [nx] , 0) , starting from 
the seed in cell (t',0), is stochastically dominated by S}_,^-,,. The estimate for 

P(il>^g(x,t)) and Borel-Cantelli imply that, with probability 1, for any rational 
(x,t) and s > 0, the event Dn^s{x,t) happens for all large enough n. 

Given now arbitrary (xo,to), pick rational (x,t) and s > so that |x — xqI < s/A 
and to < i- Then eventually c?^ A^t) < o'VAnt + ns), while by monotonicity 

^r-o](^^o)<a[;.„](nt). Weget""' 

1 

lim sup -cf'^rix,^] (nto) <i{j{x,t + s). 

n — >oo fl 
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Let t + s \to and x ^' xq, use the continuity Lemma 7.1, and conclude that now 

(7.9) holds for all (a;,t) outside a single exceptional event of zero probability. 

It remains to prove that (7.9) holds a.s. for a fixed rational {x,t). Pick rational 
< si < So, with the intention that sq \ in the end. Set Bi = (t + si) ■ Bq. 
Define 

(7.10) Cn = max {<(0) + Zf_](nt)} . 

Lemma 7.2. With probability 1, cr^l^int) = ^^ for all large enough n. 

Proof. One way to guarantee the equality crr^j-Ant) = $,n is to require that Zy.(nt) 
= — oo for all V ^ [nx] + nBi. By Lemma 4.1, distributionwise this is equivalent 
to requiring that [nx] ^ B'"{nt) for all v ^ [nx] + nSi, where B'" is a first-passage 
percolation cluster starting from a seed at site v, defined as in (2.2) in terms of Z^ . 
Switching to complements, 

< P( [nx] G B^{nt) for some v ^ [nx] + nBi ) 
= P( G B'"{nt) for some v ^ nBi ) 

< Yl PiOeB^int)) 

= J2 P{nO,v)<nt) 

< Ci exp(— C2n) 

for constants Ci = Ci{t,si), by Corollary 4.1. The conclusion now follows from 
Borel-Cantelli. D 

Let {Vi : 1 < i < ni} be a collection of closed neighborhoods whose interiors 
cover the compact set x + Bi, and such that each 

(7.11) V,Cy, + {so/2)Bo 

for any yi G Vi, each Vi lies inside x + {t + sq) ■ intSo, and satisfies assumption 
(2.6). Since the interiors cover x + Bi, we have 

(7.12) [nx] +nBi C [jnV 

i=l 

for large enough n. Pick yi G Vi such that 

(7.13) g('-^yi]<Mg(^+so. 
^ ^ \t + soJ -yev.^lt + so/ 
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By (7.11), Corollary 4.2, and Borel-Cantelli, the following holds with probability 1: 
for large enough n, 

(7.14) Z\^y^\nso) o e.r^so > for all v G nV,, 

for alH = 1, . . . , TO. In words: Start ballistic deposition processes from seeds in cells 
([n^i], 0) at Poisson process time —usq. If n is large enough, at Poisson process time 
each of these processes has grown sufficiently to cover its piece nVi. Henceforth 
assume that we are on this event of full probability, and that n is large enough for 

(7.14) to hold. 

If we define new processes Zu'''{s) = Zu {s + usq) o 6-nso-, (7.14) gives the 
inequality 

(7.15) Zl{s) < Z["^^](s) for all v G nV^ and u e Z'^ 

at time s = 0, for alH = 1, . . . , to. The monotonicity Lemma 3.2 then ensures that 
(7.15) holds at all times s > 0, and we get 

hn= max !^a:{0) + -ZrJnt) 
^ ^ n ve[nx]+nBi[n n ^"^^^ 

(7.16) 



< max <; - ■ max (j,7(0) + -Zl^^^^Unt + nso) o 6. 

l<i<m n vEnVi U ^""^J 



■nso 



Now let n ^ oo, apply assumption (2.6), Theorem 4, and the choice (7.13) of yi. 
The limit in Theorem 4 can be taken because yi & Vi (^ x + {t + sq) ■ intSo- 

1 I / X — 1J ' 

limsup -^n < max <^ sup V'o(y) + {t + so)g( -— — - 
^^■^^^ < sup L,(y) + (t + so)g(^)\+so 

yex+(t+so)-intBo L \t + Sq J ) 

= ij{x,t + so) +So. 

The argument can be repeated for arbitrarily small sq > 0. Let sq \ 0, use the 
continuity of t/j (Lemma 7.1), and then Lemma 7.2 to conclude that (7.9) holds a.s. 
The strong law of Theorem 1 under Assumption A is thereby proved. 

Under Assumption B, inequality (7.7), Lemma 7.2, and (7.16)-(7.17) give, almost 
surely, 

1 1 

M < lim inf -a?!^, (nt) < lim sup -ajl^i (nt) < ipix, t + sq) + sq 

for any M < il){x, t) and any sq > 0. This proves the statement in Theorem 2 under 
Assumption B. The argument is the same for the weak law under Assumption C. 
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